Shielding quantum discord through continuous dynamical decoupling 
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This work investigates the use of dynamical decoupling to shield quantum discord from errors introduced by 
the environment. Specifically, a two-qubits system interacting with independent baths of bosons are considered. 
The initial conditions of the system were chosen as pure and mixed states, while the dynamical decoupling is 
achieved by means of continuous fields. The effects of the temperature on the shielding of dynamical discord 
is also studied. It is found that, depending on the initial state of the system, protecting dynamical discord via 
dynamical decoupling technique may not be advantageous compared to the free-evolution of the system below 
a certain time, which depends on the frequency of the protective field. Above this time, the shielding through 
dynamical decoupling becomes effective in comparison to the system free evolution. It is also shown that the 
time for which the shielding of the quantum discord becomes effective decreases as the temperature increases. 



I. INTRODUCTION 

One of the major obstacles in building an efficient quan- 
tum computer is the presence of errors introduced by the en- 
vironment. For over twenty years, different methods have 
been developed aiming to protect the quantum information, 
such as quantum error-correction codes [1], decoherence-free 
subspaces or subsystems [2] and dynamical decoupling [3- 
6]. Among these strategies, the dynamical decoupling has the 
advantage of not requiring the use of extra qubits to encode 
the logical qubits, actively protecting quantum information 
by means of a sequence of ultrafast pulses or high frequency 
fields. 

Quantum entanglement has been considered one of the 
main concepts concerning the measurement of quantum cor- 
relations. However, several alternative measures of quantum 
correlations have drawn considerable attention in the last few 
years. Among these plethora of new correlation measure- 
ments [8], quantum discord [9, 10] has been one of the most 
employed. Quantum discord (QD) is intimately connected to 
the entanglement of formation [11-13], the conditional en- 
tropy [13] and the mutual information. Moreover, quantum 
discord is related to many important protocols as the distri- 
bution of entanglement [14], quantum locking [15], entangle- 
ment irreversibility [16], and many others [17]. 

For some specific conditions [18], discord has the peculiar- 
ity of not being affected by the environment during an initial 
finite time. Such a finite time is known as transition time be- 
cause it delimits a sudden transition from classical to quantum 
decoherence regime. Recently, Luo et al. [19] investigated 
the possibility of extending the transition time by using bang- 
bang pulses to control the decoherence of two qubits coupled 
to independent spin baths. Their results demonstrated that 
the transition time can indeed be extended through protection 
techniques. In this work, we theoretically investigate the pro- 
tection of the QD by means of the dynamical decoupling (DD) 
technique. We consider a situation where the DD is achieved 
by using continuous fields, which is a more realistic situation 
than using instantaneous pulses [20]. The interest of the scien- 



tific community has been increasing for this kind of protection 
and very recently it has been experimentally implemented to 
protect the quantum information in nitrogen vacancy systems 
[21, 22]. It has been shown that using this kind of protection, 
the coherence time can be prolonged by about twenty times 
[22]. Furthermore, it can be applied to a wide variety of phys- 
ical systems including trapped atoms and ions, quantum dots 
and nitrogen-vacancy centers in diamond. 

In order to probe the protection method, we solve the Red- 
field master equation taking into account the interaction of 
each qubit with their own boson bath. We study the effects 
caused by temperature and two kinds of initial conditions: 
pure and mixed states. We show that protecting QD by means 
of DD technique can be tricky depending on the initial con- 
dition and the desired time scale. Contrarily to entanglement 
and fidelity, quantum discord cannot be protected by means 
of DD for some particular situations. This result shows a pe- 
culiar aspect concerning the protection of the quantum dis- 
cord: sometime keeping the free-evolution of the qubits may 
be advantageous in comparison to the protection by DD below 
a certain specific time. Such a result is very different from 
the one obtained through instantaneous pulses approach [19]. 
This difference occurs because instantaneous pulses applied to 
the system do not break the symmetry of the reduced density 
matrix, which is an essential ingredient to keep QD intact be- 
fore the transition time [18]. On the other hand, a continuous 
field applied to the system induces a dynamics that breaks the 
symmetry of the reduced density matrix, thereby destroying 
the sudden transition. 



D. DYNAMICAL DECOUPLING 

The Hamiltonian for a two-qubit system and the environ- 
ment can be written as 

H = H + H E + H int , (1) 

where Hq is the environment-free Hamiltonian, which in- 
cludes the qubits Hamiltonian and the time-dependent fields 
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applied to protect the quantum discord. H E is the Hamil- 
tonian describing the environment and the Hamiltonian H int 
takes into account the interaction between the qubits and the 
environment. Here, 
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represents two independent baths of harmonic oscillators (one 
for each qubit) with uj l k given by the frequency of the normal 

mode k of the thermal bath and a\ (ajj) is the usual annihi- 

lator (creator) operator for the i-th qubit,/.e. at (a , k ) is the 
operator that annihilates (creates) a bath quantum in mode k 
of the i-th qubit. 

The interaction Hamiltonian H nt describing two qubits 
coupled to their baths is written as 
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where £j, for i — 1,2, are operators that act on the environ- 
mental Hilbert space, ay are the Pauli z matrix acting on 
i-th qubit . This model is called a dephasing noise and it is 
a typical class of errors in quantum dots [23, 24] and Joseph- 
son junction [24]. Explicitly, d = Bi + B[, Bi = ^2 k g l k a\, 
where gj. is a complex coupling constant for fc-th normal mode 
with frequency dimension. Henceforth, we use units of ft = 1. 

For the sake of simplicity, we assume there exist a global, 
static control field that exactly cancels the intrinsic system 
Hamiltonian, i.e. Hq represents only the active protection and 
it is given by a continuous applied field defined as 
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where cj = 2tt /t c , t c is the period of Uo(t) = exp(-iHot), 
and n x is an integer number. The general prescription of the 
dynamical decoupling [3, 7] is based on finding an controlled 
unitary operation such that 



uUt)H int U (t)dt = 0. 



(5) 



This is a sufficient condition to shield the effect of the envi- 
ronment over the system and, as exposed in ref. [25], perfect 
shielding can be reached when n x tends to infinity. 



III. MASTER EQUATION 

In this paper we assume that the interaction between the 
qubit and its environment is sufficiently weak such that linear- 
response theory holds. Thus, in the interaction picture, the 
master equation is written as 
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dt 



= - f dt'Tr E {[H r (t), [#/(*'), M0)p/(*)]]} , (6) 



where Hj(t) reads 



ff/(t) = U%(t)U E (t)H int U E (t)Uo(t), 



(7) 



and 



Uo(t) = exp -i 



.Hot 



(8) 



(9) 



Pe(0) is the initial density matrix of the thermal bath, which 
is chosen as the Feynman- Vernon state, 



PE (0) = ^exp{-j3H E ), 
where Z is the partition function given by 

Z = Tr E [exp(-0H E )] , 



(10) 



(11) 



/3 = l/(fcsT), ks is the Boltzmann constant, and T is the 
temperature of the environment. Defining U E {t)CiU E (t) = 

Ci(t) and Ai(t) = Ul(t)a ( z l) U (t), we can write the Hj 
Hamiltonian in the interaction picture as follows, 



H J (t)=A 1 (t)£ 1 (i) + A 2 (t)£ 2 (i). 



(12) 



Thus, substituting Eq. (10) and Eq. (12) into the master equa- 
tion Eq. (6), we obtain 
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where T>i(t, t') = T l (t- t') + %(t - if) with 

Xi(t-t') = Tr E {Bi(t)p E (t)Bltf)} 

= J2\9i\ 2 nieM-^Ut-t')} (14) 



and 



Ti(t-t') = Tr E {Bl(t)p E (t)Bi{t')} 

= E ISfcl'K + 1) exp[-i4(f - t% (15) 



with Bi(t) = U E (t)BiU E (t), p E (t) = U E (t)p E (0)U E (t), 
and n\ is the average occupation number of mode k of the i-th 
qubit, n\ = l/[exp(/3a;ji. — 1)], Since each qubit is subjected 
to identical environments, the i index of g k , n k and uj k can 
be omitted. Thus, defining the spectral function as J{oS) = 
Ylk \9k\ 2 S(uj — Wfc), we can replace the above summations by 
the following integrals 



and 



X(t) = / du> J(uj)n(uj) exp(— iwt), (16) 



T(t)= aW(w)exp(iwi)[n(w) + 1]. (17) 
Jo 
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To illustrate our results, we assume that the environment has 
an ohmic spectral density, J(uj) = tjuj exp(—u/u) c ), where oj c 
is a cut-off frequency and r} is a damping constant. Therefore, 
we can explicitly evaluate the above integrals which yield 



V{t,t') = 
2// 



[1 + iu c (t - t')} 2 



+ ^Re |V lst )(l + l/GSwc) - i(t - f )/£)}, (18) 

where »3/( lst ) is the first polygamma function. Here, Eq. (18) 
accounts for both the interaction with the electromagnetic 
vacuum and with photons that could be present in the envi- 
ronment. Indeed, the first term of the righthand side repre- 
sents the vacuum and the second term accounts for the ef- 
fects of finite temperature. In particular, for T = we have 
V(t,t') = r)cjl/[l + iu) c {t-t')} 2 . 



IV. FIDELITY, CONCURRENCE, AND QUANTUM 
DISCORD 



In this section we present the definition of three distinct 
measures for quantum systems. 



C. Quantum Discord 

The quantum discord, was independently defined by Han- 
derson and Vedral [9], and Ollivier and Zurek [10] about ten 
years ago. It was proposed as a measure of the quantum cor- 
relations and, given two subsystem A and B, it is defined as 
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where Iab = Sa + Sb — Sab is the mutual information, with 
Sa = — T?(pa log Pa) given by the von-Neumann entropy 
of the subsystem A and similarly to B and AB. J AB is the 
classical correlation, which can be explicitly written as 



Jt~ n = max 
AB {n fc } 



Sa - ^2pkS{p A \k) 



(22) 



where p A \k = Tr B {I\. k p A B^k) /Tr ab {likPAB^k) is the re- 
duced state of A after obtaining the outcome k in B. Here, 
{n^} is a complete set of positive operator valued measure- 
ment that results in the outcome k with probability p k . The 
quantum discord measures the amount of the mutual informa- 
tion that is not locally accessible [13, 28] and generally is not 
symmetric, i.e. 8^ B ^ 8 BA . 



A. Superfldelity 

To calculate the dissipative dynamics fidelity, we use an al- 
ternative measure known in the literature as superfldelity [26]. 
We calculate the superfldelity between two distinct (possi- 
bly mixed) density matrix state, p(t) and a(t). Here, p{i) 
represents the open quantum system dynamics, where the 
two qubits interact with the environment and a(t) gives the 
ideal dynamics, whose interaction with the environment is dis- 
able. The superfldelity is a function of linear entropy and the 
Hilbert-Schmidt inner product between the given states. Re- 
markably, the superfldelity is jointly concave and satisfies all 
Jozsa's axioms. Explicitly it is given by: 

F(p, a) = Tv(pa) + y/l - Tv(p 2 )^/l - Tv(a 2 ), (19) 
and it is valid for pure or mixed states. 



B. Concurrence 

To measure the entanglement between the two qubits, we 
use the well know concurrence [27]. It is defined as the maxi- 
mum between zero and A(t), 

A(t) = Ai - A 2 - A 3 - A 4 , (20) 

where Ai > A2 > A 3 > A4 are the square roots of the eigen- 
values of the matrix p(t)cr y ® a y p*(t)cr y (g) a y , where p*(t) is 
the complex conjugation of p(t) and <j y is the y Pauli matrix. 



V. DISSIPATIVE DYNAMICS 

In order to illustrate the singular feature involving the pro- 
tection of the quantum discord, we study two kinds of initial 
conditions, pure states and the mixed ones. Additionally, we 
consider the effects of finite temperatures on the quantum dis- 
cord. Moreover, the quantum discord dynamics is compared 
to the concurrence and the fidelity of the quantum systems. In 
the following calculations, we set the parameters of Eq. (18) 
as being 77 = 1/16 and u c = 2tt. 



A. Pure-state initial condition 

In order to analyze the dissipative dynamics for pure states, 
we choose for the initial condition a maximal quantum corre- 
lated state PAB {0) = -^(|00)(00| + |11}(11|). 

Figure (1) shows the concurrence (a) and the fidelity (b) 
dynamics when the control fields are turned off (n x = 0) 
and on (n x = 4,6,8) for a zero temperature environment, 
T = OK. The solid line corresponds to n x = in Eq. (4) 
and the dashed, dotted, and dot-dashed lines give the concur- 
rence and the fidelity when the control fields are turned on 
with n x = 4,6,8, respectively. When the frequency of the 
protective field increases both the fidelity and the concurrence 
increase [29]. 

In Figure (2), quantum discord is plotted considering the 
same parameters used for fidelity and entanglement. The solid 
line corresponds to n x = in Eq. (4) and the dashed, dotted, 
and dot-dashed lines give the quantum discord when the con- 
trol fields are turned on with n x = 4, 6, 8, respectively. For 
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FIG. 1: (Color Online) Concurrence and fidelity as a function of time 
considering a pure initial state and T = OK. (a) The solid (black) 
line gives the concurrence when the control field are turned off and 
the dotted (red), double dot-dashed (blue), and dot-dashed (green) 
lines give the concurrence when the control fields are turned on, us- 
ing n x = 4, 6, 8 respectively, (b) The solid (black) line gives the 
fidelity when the control field are turned off and the dotted (red), dou- 
ble dot-dashed (blue), and dot-dashed (green) lines give the fidelity 
when the control fields are turned on, using n x = 4, 6, 8 respectively. 
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FIG. 2: (Color Online) Quantum Discord as a function of time con- 
sidering a pure initial state and T — OK. The solid (black) line gives 
the quantum discord when the control field are turned off and the 
dotted (red), double dot-dashed (blue), and dot-dashed (green) lines 
give the quantum discord when the control fields are turned on, using 
Tlx = 4, 6, 8 respectively. 



the initial pure state case, the concurrence, fidelity, and quan- 
tum discord presents a similar behavior, i.e. as the frequency 
of the protective field increases, the protection of the quan- 
tum information is enhanced [29]. It is important to empha- 
size that, despite the restrictive initial state, we point out that 
the fidelity, entanglement, and quantum discord protection is 
reached for all pure initial state and, qualitatively, in the same 
way, i.e. the higher the frequency of the protective field, the 
better the protection. As shown in the following, the quantum 



discord does not necessarily behaves in the same qualitative 
way as fidelity and concurrence with respect to the applica- 
tion of the protective field for mixed states. 

B. Mixed-state initial condition 

For mixed states, we choose as initial condition a maximal 
mixed marginal state. The initial condition for this situation 
is given by pab(0) = 0.8|$}{$| +0.2|*)(*|, where |$) = 
^(|00) + |ll))and|*) = ^(|01) + |10)). 

In Fig. (3), we show the concurrence and the fidelity dy- 
namics for some frequencies of the protective field. There is 
no qualitative difference involving the protection of the con- 
currence and the fidelity and the greater the frequency of the 
protective field, the better the protection. On the other hand, 
we observe a very particular behavior when dealing with the 
quantum discord. It can be noticed in Fig. (4) an abrupt de- 
crease, the sudden transition, in the quantum discord around 
t = 0.4t for the non-protected case. Up to this sudden tran- 
sition, the results obtained with the protection turned on are 
less efficient in maintaining the quantum discord than with 
the protection turned off. Actually, for t < 0.4r, the protec- 
tive dynamics is equivalent to the non-protected one, only in 
the limit that the frequency of the protective field goes to in- 
finity. However, there is a finite time when the protection of 
the quantum discord surpass that of the non-protected case, 
given by the crossing of the dotted lines with the solid line in 
Fig. (4). The effectiveness of the protective field only occurs 
above these times t e (n x ), which are greater than the time of 
sudden transition of the quantum discord t e {n x ) > 0.4r. 

This peculiar aspect involving the protection of the quan- 
tum discord occurs because the plateau given in the free evo- 
lution exist for some special symmetries of the initial state 
[18, 30]. The dynamical decoupling strategy is an active 
method of protection and needs to continuously modify the 
density matrix during the dynamics. On the other hand, the 
density matrix reaches its original and special symmetry just 
at some periodic times. This necessary density matrix modi- 
fications involving the protection by means of dynamical de- 
coupling decreases the quantum discord and reaches a value 
less than the initial one. However, when the quantum discord 
dynamics reaches the sudden transition time, it starts to de- 
crease while the protective scheme becomes more relevant. In 
the following, we show how a finite temperature environment 
influences this scenario. 

The first relevant observation involving finite temperature 
environments is that the interaction with the environment be- 
comes more deleterious. The density matrix coherence decays 
faster and the same occurs to the entanglement, the fidelity and 
the quantum discord. As a consequence of this faster decay, 
the quantum discord plateau diminish. Furthermore, the time 
for effectiveness of the protective field becomes smaller as the 
temperature increase. In other words, the effective protection 
occurs for smaller time scales, as the temperature increases. 

In Fig. (5), we study the temperature effects on the protec- 
tion method. We present the time for effectiveness of the pro- 
tective field for different temperatures, as a function of the fre- 
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FIG. 3: (Color Online) Concurrence and Fidelity as a function of 
time considering a mixed initial state and T = OA", (a) The solid 
(black) line gives the concurrence when the control field are turned 
off and the dotted (red), double dot-dashed (blue), and dot-dashed 
(green) lines give the concurrence when the control fields are turned 
on, using n x = 4, 6, 8 respectively, (b) The solid (black) line gives 
the fidelity when the control field are turned off and the dotted (red), 
double dot-dashed (blue), and dot-dashed (green) lines give the fi- 
delity when the control fields are turned on, using n x — 4, 6, 8 re- 
spectively. 
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fore, for higher temperatures, the protective scheme becomes 
efficient sooner. This characteristic occurs because when the 
temperature increases, the density matrix coherence decays 
faster and, consequently, the time of sudden transition is re- 
duced. Thus, the plateau of quantum discord is smaller and 
the effectiveness of the protective scheme is faster. 
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FIG. 5: (Color Online) The time for efficiency as a function of 
the protective field frequency for different temperatures (T = 
OK, 0.1K, 0.2K). 
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FIG. 4: (Color Online) Quantum Discord as a function of time con- 
sidering a mixed initial state and T — OK. The solid (black) line 
gives the quantum discord when the control field are turned off and 
the close points (magenta), dotted (red), double dot-dashed (blue), 
and dot-dashed (green) lines give the quantum discord when the con- 
trol fields are turned on, using n x = 2, 4, 6, 8 respectively. In the 
inset we plot the time for efficiency as a function of n x . 

quency of the protective scheme. As one can observe, t e (nx) 
becomes smaller when the temperature is increased. There- 



VI. SUMMARY 

In summary, we have shown that even that the DD pro- 
tective scheme increase the fidelity or the entanglement, the 
quantum discord can decrease when compared to the evolu- 
tion without protection. Indeed, the protection of the quantum 
discord can be a very difficult task depending on the initial 
condition and the time scale of interest. It was shown that 
this peculiar aspect involves the time scale desired to protect 
quantum information and that, for a time less than the time 
of sudden transition, the protection is not efficient. We have 
compared two kinds of initial condition, the pure one and the 
mixed one and included finite temperature effects. We have 
verified that the greater the temperature, the faster the effi- 
ciency of the protective method is achieved, which evidences 
the the relevance of the temperature in the protection of QD. 
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